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Theoretical calculations, computer simulations and experiments indicate the possible existence 
of a ferromagnetic liquid state, although definitive experimental evidence is lacking. Should such 
a state exist, demagnetization effects would force a nontrivial magnetization texture. Since liquid 
droplets are deformable, the droplet shape is coupled with the magnetization texture. In a thin- 
film geometry in zero applied field, the droplet has a circular shape and a rotating magnetization 
texture with a point vortex at the center. We calculate the elongation and magnetization texture 
of such ferromagnetic thin film liquid droplet confined between two parallel plates under a weak 
applied magnetic field. The vortex stretches into a domain wall and exchange forces break the 
reflection symmetry. This behavior contrasts qualitatively and quantitatively with the elongation 
of paramagnetic thin films. 



I. INTRODUCTION 



The study of a possible ferromagnetic liquid state is a problem of considerable interest. In such a spontaneously 
magnetized liquid state, long range magnetic order would exist in the liquid without application of any external field. 
The existence of such a liquid state has been indicated by mean field calculations (l]-^| and computer simulations on 
strongly polar fluids 

Experiments to observe ferromagnetism in liquids with strong magnetic interactions such as ferrofluids [fTo|| have 
failed so far because the liquids freeze Jll],[l2| or phase separate |l3| well above thepredicted low temperatures for 
the onset of spontaneous magnetization. Experiments on super-cooled Co-Pd alloys B] do indicate the possibility of 
ferromagnetism in liquids. In this case it is the strong exchange interaction and not the dipole interaction that would 
cause the spontaneous magnetization. The experimental evidence regarding Co-Pd is still not conclusive. 

Although the existence of a ferromagnetic liquid state is yet to be confirmed experimentally, spontaneous polariza- 
tion coupled with other order parameters has already been observed. Some electrically polarized liquid crystals |T^ 1 
show a helical ordering of the dipole moments in the liquid. In superfluid 3 He the magnetic moment couples to the 
superconducting order parameter Jl6| . Many superfluid 3 He phases are therefore also magnetically ordered. 

It is interesting to consider the magnetization texture (spatial variation of the orientation of magnetization) inside a 
droplet of such a ferromagnetic liquid |l7j . The magnetization texture likes to avoid poles to minimize its energy. 
However, this leads to defects inside the texture. For example, a rotating magnetization texture with cylindrical 
symmetry inside a sphere avoids all poles but has a vortex line running through the center. Near the vortex of 
such a texture the magnetization is topologically unstable [ [l9| and might escape into the third dimension pC| ] with 
a nonzero component along the vortex line. Whether this happens depends on the balance between demagnetizing 
and vortex energies. Simulated annealing of the magnetization inside a cubic box suggests that replacing vortices 
with point defects may be favorable Q . Any defect is likely to have a system-shape dependent energy cost causing a 
deformable liquid droplet to deviate from a spherical shape. The complete calculation of the shape of an unconfined 
ferromagnetic liquid droplet in three-dimensions coupled with the calculation of its magnetization texture remains an 
unsolved problem. 

This problem has a simple solution in two-dimensions in zero field. The magnetization texture inside any soft (zero 
anisotropy) ferromagnetic solid thin film is given by van den Berg's algorithm which avoids all poles, and thus 
all magnetostatic energy, at the expense of a domain wall through the film. A liquid droplet, which can change its 
shape, prefers a circular shape to minimize its surface energy. The magnetization lines inside a circle form concentric 
circles according to van den Berg's algorithm. For a circular shape the domain wall energy is also minimized because 
the domain wall shrinks to a point vortex. The circle thus solves the coupled texture and shape problem in zero field. 

Our goal is to analyze the shape of a ferromagnetic liquid droplet under weak applied fields. Since the texture and 
shape of a three-dimensional drop is not precisely known at present, we concentrate on the two-dimensional case of 
a droplet confined between two parallel plates with spacing much smaller than the droplet diameter. When a field is 
applied parallel to the plates the magnetization texture distorts to exclude the magnetic field from the bulk of the 
droplet. Bryant and Suhl |2|] calculated the texture of solid circular and elliptic thin films under applied field. We 
adapt their result to ferromagnetic liquids by letting the shape vary as the field is applied. We find that the droplet 
elongates to reduce its magnetostatic energy. The equilibrium shape is reached when the magnetostatic energy plus 
the surface energy is minimized. 
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In Sec. || we calculate elongation of a ferromagnetic th in film as a function of its undeformed radius, thickness, 



the saturation magnetization and the applied field. In Sec. Ill A we contrast this behavior with that of paramagnetic 
thin films under applied fields and discuss why this contrasting behavior occurs. In Sec. Ill B we discuss the sym- 
metry breaking of the droplet shape under exchange forces and its dependence on the applied field. Finally, Sec. I\ 
summarizes our results. 

II. CALCULATION 

Consider a homogeneous model liquid such that every point in the liquid has a saturation magnetization M s with 
no anisotropy. Let a droplet of such a liquid in a spontaneously magnetized state be confined between two parallel 
plates as shown in Fig. 1(a). The thickness spacing A between the plates is kept much smaller than the diameter 
of the undeformed droplet so that the droplet has a very small aspect ratio p = A/2ro. Figure 1(b) shows a droplet 
elongating under applied field. To a good approximation the droplet shape is an ellipse for small elongation. We 
define the elongation as 

e = a/6-1, (1) 

where a and b are the semimajor and semiminor axes of the ellipse. We choose a coordinate system with the field in 
the x direction and z normal to the plates. The origin is at the center of the droplet. 

Under zero external field this thin film droplet has a circular shape. The magnetization texture M(r), given by van 
den Berg's algorithm J2l| , forms circles concentric with the droplet boundary. Because this texture is divergence-free 
and everywhere tangent to the surface, there are no magnetic poles Jl8[ and hence no demagnetizing field. This 
texture therefore achieves the lowest possible magnetostatic energy. When an external field is applied, the circular 
magnetization texture distorts to reach a new equilibrium configuration. The new texture exhibits poles, making the 
magnetostatic energy shape dependent. The droplet elongates to reduce this shape dependent magnetostatic energy 
and an equilibrium shape is reached when the magnetostatic plus surface energy is minimized. We n eglect magnetic 



exchange energy because it is small relative to the magnetostatic energy as discussed in Sec. Ill B . Our goal is to 
calculate this elongation as a function of the applied field and the various parameters for the droplet. 

To calculate the magnetostatic energy of the droplet we make use of the calculations by Bryant and Suhl |^2| for 
the magnetization texture of elliptic solid ferromagnetic thin films under applied field. Their calculations for the 
magnetization texture are for fixed shape. We let the shape of the thin film droplet vary and calculate its energy as 
a function of the shape for small elongation. The elongation is then calculated by minimizing the total energy with 
respect to elongation. 

Bryant and Suhl calculate the equilibrium magnetization texture of an elliptical ferromagnetic thin film in applied 
field by drawing analogy to the case of conductors in applied electric field. A ferromagnet is equivalent to a conductor 
because it has effectively an infinite permeability. When a conducting film is placed in an applied electric field, it 
expels the electric field by creating an induced charge density on its surface. A soft ferromagnet will similarly expel 
magnetic field if a magnetization texture can be found exhibiting the necessary surface pole density. These charges 
&m(%, y) arise from the components of M normal to the top and bottom surfaces of the droplet and there are no poles 
within the bulk of the droplet. However, for a thin film we treat the charge as if it permeates through the film and 
the magnetization texture is restricted to the x — y plane so that M z = 0. The magnetization texture then satisfies 
the following pseudo-divergence equation: 

„ , , dM x dAL , w * 

W-M = ^ + -^ = -p M = -a M (x,y)/A. (2) 

The above equation is exact in the limit of zero thickness. It remains a good approximation provided the the film has 
a small aspect ratio (A/2ro << 1), because strong demagnetizing fields prevent the magnetization from tilting far 
out of the x — y plane. The partial differential Eq. can be integrated, in principle, to find a solution for M(r) if it 
exists. 

For an elliptic thin film, the charge density that expels the field can be written in a convenient form as 

VM{x,y) = — H , (3) 

na 2 (l-(x/a) 2 -(y/b) 2 ) 1 

where 



E x = Ho \ = =r- (4) 

\8M s bA(K-E)J 
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is the reduced applied field. K and E are complete elliptic integrals of the first and second kind, respectively, of the 
argument (1 — 6 2 /a 2 ) 1 / 2 . For E x < 1, Eq. (0) can be integrated to find the magnetization texture which expels the 
magnetic field. A C-shaped domain wall (see Fig. 2) appears inside the droplet for nonzero field. For E x > f , the 
domain wall intersects the boundary of the film and the field penetrates the interior of the film. 

We use this result to calculate the magnetostatic energy of a liquid droplet. Since the above calculations are for 
zero anisotropy, they also apply to liquids. We begin with the general expression for the total magnetostatic energy 
of any magnetization distribution under an applied field 



dr H • M — - / <1t II • M. (•-.) 



The first term on the right hand side of (|J) is the energy of the applied field Ho acting on magnetization M. The second 
term is the self-energy (hence the factor of 1/2) due to the magnetization interacting with its own demagnetizing field 
Hd- For E x < 1, Hd — —Ho because the magnetization expels the field, therefore the magnetostatic energy 

E M = -l f drUo-M. (6) 
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Using pm = — V • M, writing Ho = — V0o 5 an d integrating by parts, the magnetostatic energy in (||) can be rewritten 

Em = 7; J dr fa p M , (7) 

where (f>o = — Hqx is the potential due to Ho. Using pm = &m/A and Eqs. (|^) and (^) and calculating the integral 
in ^ gives 

32 o n A 9 .9X — E 1 , a 2 - b 2 
E M = ^ElM^ 2 ab 2 — = -^ a„, (8) 

Assume that the applied field is small so that e << 1. Expanding the result in (||) for small e gives 

9H 2 r 3 3 

E M = - 2 -^ { l + \e), (9) 

with corrections of higher order in e and higher order in the aspect ratio A/2ro. As expected, the leading correction 
to the energy is negative, so the droplet elongates to reduce its magnetostatic energy. Interestingly, the magnetostatic 
energy is independent of the saturation magnetization M s . This is so because, as long as E x < 1, the charge density 
Pm does not depend on M s . The factor of M s in the numerator of (tm( x > y) m (||) cancels against the factor of M s 
in the denominator of E x in Eq. (|]). Physically, the charge density of a given distribution is proportional to M s , but 
the distortion that creates pM varies inversely with M s . 

The magnetostatic energy is also independent of A to the lowest order in the aspect ratio A/2r . To understand 
this, note that the expression for the reduced applied field E x has a factor A in the denominator implying that thicker 
films can expel higher values of applied field for a given distortion. Thus a thicker film has a smaller distortion in 
its texture, for a particular value of applied field, than a thinner film. This is confirmed by the factor of the A 
in the denominator of the charge density pm- When the energy density is integrated over the volume it yields a 
magnetostatic energy independent of A. 

As the droplet elongates, its perimeter increases and is given by 

S = 27rr (l + ^e 2 ), (10) 

with higher order corrections in e. The variation of the perimeter is quadratic in e as expected because the perimeter 
should increase regardless of the sign of e. Since the area of the droplet in contact with the plates remains constant, 
the relevant surface energy is the surface energy along the perimeter, 

E s = aSA, (11) 

where a is the surface tension along the perimeter of the droplet. Note that we consider the case of 90° contact angle 
of the magnetic fluid- non magnetic fluid interface with the boundary plates. Discussion of surface energy for other 
contact angles is in Ref. |23| ]. 

Minimizing the total energy Em + Es with respect to e gives the elongation of a ferromagnetic droplet, 
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Eferro = 2 T ■ (12) 

The elongation is quadratic in the applied field since reversing the field direction should not affect the result. At 
fixed thickness A, droplets with larger radius r will elongate more. This is because magnetostatic energy in ([)]) is 
proportional to rg, whereas the surface energy is proportional to the area of the curved surface (of order ?"oA). The 
elongation, however, is inversely proportional to A. This is because the droplets with larger thickness find it easier to 
expel the applied field. The same applied field causes a relatively smaller distortion in the magnetization texture of 
a thicker droplet, and a smaller distortion of texture causes a smaller elongation. Figure 3 shows the magnetization 
texture of a droplet as it elongates under applied field. 

Determining the applied field at which the domain wall touches the boundary and field penetration occurs is 
complicated by the dependence of the reduced field E x , defined in Eq. (Q), on elongation. Basically, a given applied 
field Ho causes a greater distortion of the magnetization texture in a droplet that is elongated in the field direction 
than it would have caused in a circular droplet. This is because the demagnetizing field becomes weak in elongated 
droplets (that is why the droplet elongates, after all) so a larger charge density is needed to achieve a demagnetizing 
field H/j = H . Correspondingly, a lower applied field is needed to achieve penetration for an elongated droplet 
than for a circular droplet. 

The reduced applied field in (|4|) can be written to first order in e as 

with corrections of order e 3 . Substituting the expression for elongation ([l2|) into the above gives 

^ = -^(l + ^_). (14) 

For any finite <r, E x now depends nonlinearly on Hq. E x is still a useful quantity as it measures the distortion in the 
magnetization texture of the droplet, but its relation to applied field is nontrivial. It is therefore convenient to define 
a new dimensionless field, 

Hor ° (15) 



2ttAM s ' 

that equals E x when the droplet is circular. We also define a dimensionless surface tension 

(16) 



so that the reduced field can be written as 



2M?A' 



h* 2 



E x = h*(l + —). (17) 



Our calculations hold in the limit of small droplet distortion, h* 2 /a* << 1, and non-penetration of the field, E x < 1. 
Consider the applied field at which field penetration occurs. Solving Eq. ([jj]) for the critical applied field h* for which 
E x = 1 yields h* w 1 — 1/ct* in the limit of nearly undeformable droplets, a* >> 1. 

Figure || illustrates the simultaneous evolution of droplet shape and texture. Three stages are shown, at E x = 0, 0.5 
and 1. Parameters are chosen so that a* — I. As expected, in contrast to the values of E x , the actual applied field 
values h* become more closely spaced as the elongation grows, and h* c = 0.68 < 1. 

III. DISCUSSION 

A. Contrast with paramagnetic thin films 

Expression (|l^) for elongation of ferromagnetic droplets contrasts with the elongation of paramagnetic droplets in 
the same thin film geometry. For a paramagnetic thin film with susceptibility \i the elongation has been previously 
calculated |E3t|, 



4 



X 2 HjA Cr 

Cpara = m —7- , (18) 

a A 

where C is a constant. 

To understand this contrasting behavior we provide the following simple argument on dimensional grounds. The 
elongation is a ratio of the magnetostatic and surface energies. Let Q be the net charge on the portion of the droplet 
with x > 0. For x < the net charge is —Q. In both the ferromagnetic and paramagnetic cases, the charges on 
average are separated by distances of order r , and therefore the magnetostatic energy goes as Q 2 /r . Dividing this 
by the surface energy (of order tqA<j) gives the elongation 

<}1 (19) 



❖A' 



where "~" indicates proportionality. In the case of a ferromagnetic droplet Q ~ H^r^ because the charge density pM 
is distributed in the bulk of the droplet (volume ~ A) but is inversely proportional to A. This explains why the 
elongation of a ferromagnetic droplet increases rapidly with r and decreases with A in Eq. (|l2|). In the case of a 
paramagnetic droplet, the charges are distributed mainly along the curved surface of the droplet, so the total charge 
Q is therefore proportional to xHqTqA. Our simple dimensional argument gives e ~ x 2 HqA/<j, which reproduces 



Eq. (18) except for the dimensionless ln(r /A), and explains why the elongation in the paramagnetic case depends 
only weakly on ro and increases strongly with A. 

Although theory and computer simulations predict the existence of a ferromagnetic liquid state there is still not 
much experimental evidence for it. The contrasting behavior of ferromagnetic thin films and paramagnetic thin films 
under applied fields could serve as a useful test to detect such a state. This contrasting behavior can be neatly 
captured in a picture if one plots e/A vs the inverse aspect ratio, 2ro/A at a fixed field value. Figure 4 shows such 
a plot for experiments on ferromagnetic and paramagnetic droplets of a hypothetical fluid at Hq = 20 G. The fluid 
has a susceptibility x — 5 m the paramagnetic phase, surface tension a = 60 dynes/cm and M s = 500 G, typical of 
ferrofluids. We take the droplet dimensions 2r = 120 /jm and A = 1.2 /im so that p — 2r /A = 10~ 2 and a* = 1. 
The value Ho = 20 G corresponds to h* — 0.16 for the droplet in the ferromagnetic phase. Recall that for a* = 1 
field penetration occurs at h* — 0.68. 



B. Symmetry breaking of shape under exchange interaction 

Despite the asymmetric magnetization texture, the droplet shape in Fig. 3 has reflection symmetry about the 
field direction. Interactions of the magnetic poles with the applied field and with each other create the forces that 
deform the droplet. Since, for a symmetric shape the induced charge density is symmetric, the droplet deformation 
is symmetric. The deformed magnetization texture is asymmetric because the direction of rotation of magnetization 
lines (counterclockwise in Fig. 2 and 3) breaks the left-right symmetry. If the exchange interaction due to spatial 
variation of magnetization were to be included, the droplet shape would become asymmetric to lower its exchange 
energy. 

For an isotropic medium, the exchange energy density can be written as p3 

Uex = IJ^dM^ (20) 
2 dxi dxi ' 

where a is the exchange constant and the summation convention is employed. The magnetization texture for a droplet 
under applied field breaks the left-right symmetry (see Figs. 2 and 3). The breaking of the symmetry is most evident 
at the domain wall where Uex is singular. The exchange interaction creates forces that break the left-right symmetry 
of the droplet shape. 

The main contribution to the exchange energy of the droplet comes from near the domain wall, where the magne- 
tization has an apparent discontinuity of order M s . The domain wall actually spreads to a finite width w to lower 
its exchange energy at the cost of acquiring a demagnetizing energy. Since the magnetization changes by an amount 
of order M s over a distance w, the domain wall exchange energy density is roughly aM 2 /w 2 . For a Bloch domain 
wall |^5| with w << A the demagnetizing energy density is of order M 2 w/A. Minimizing the total energy per unit 
domain wall length with respect to the width gives w of order (aA) 1 / 3 and energy density of order a 1 / 3 M 2 /A 2//3 . 

For E x << 1, symmetry arguments and numerical calculations suggest that domain wall arc-length varies as E 2 r . 
The total energy of the domain wall can be estimated by multiplying the domain wall energy density with the domain 
wall arc-length (E 2 ro) and cross section (uA) to get 
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E wall ~ a^MlEl^r, = a^H^. (21) 

The ratio of the domain wall energy to the magnetostatic energy ([)]) goes like (a/A 2 ) 2 / 3 . Therefore, the relative 
importance of exchange effects diminishes for large thickness A. 

To understand how exchange interaction breaks the symmetry of the shape, we analyze the exchange forces on a 
circular disk in an applied field (Fig. 2). We write the shape of the droplet as a perturbed circle breaking the left-right 
symmetry, 

ip 2 

r(6) =ro(l--y + ^ sin 36»). (22) 

Here 9 is the angle measured counterclockwise from the direction of the applied field in the plane of the droplet. The 
sin 38 term is the lowest harmonic in 9 which breaks the left-right symmetry of the shape. The perturbation keeps 
the area constant to ensure that the volume of the droplet does not change. Assuming analytic variation with respect 
to the perturbation and the applied field, the domain wall energy takes the following form for small ip and E x , 

E waU = a 2 / 3 M 2 E 2 x r A 2 / 3 (A + Bif,E x ), (23) 

where A, B are dimensionless constants. The order ip term is multiplied by E x because under simultaneous reversals 
of ip and E x the domain wall energy remains the same. 

The surface energy of the droplet varies quadratically with ip, 

E s = 2nr Q Aa(l + 2iP 2 ). (24) 

The magnetostatic energy (|^) of the droplet has an order ip 2 correction that can be neglected because it is higher 
order in the applied field. Minimizing Es + E wa ii with respect to ip gives 

Ba 2 / 3 M 2 E 3 rn 

* = - 8,AVV • (25) 

We conclude that for weak applied fields the symmetry breaking in shape will be much smaller than the elongation 
of the droplet which goes like E 2 . The direction of symmetry breaking (sign of ip) remains undetermined, since it 
depends on the sign of B. We expect the direction of symmetry breaking to shorten the length of the domain wall. 



IV. CONCLUSION 



We calculate the magnetization texture and elongation of a ferromagnetic thin film under weak applied fields. The 
point vortex favored at zero field elongates into a curved domain wall in applied field. We find that the droplet 
elongation is proportional to the square of the radius of the undeformed droplet and inversely proportional to its 
thickness. When domain wall energies are considered, the ferromagnetic droplet may break reflection symmetry. 
These results contrast with that of paramagnetic thin films. We propose this different behavior under applied field as 
one possible way to detect the ferromagnetic state in a liquid. 
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FIG. 1. (a) A side view of a ferromagnetic liquid droplet confined between two parallel plates, (b) A top view of the droplet 
elongating under applied field. The dashed line shows the undeformed droplet. 
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E x = 




FIG. 2. The magnetization texture of a solid circular ferromagnetic thin film under applied field [22]. E x is the reduced 
applied field defined in Eq. (W). The magnetization texture in zero field is going counterclockwise. 
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FIG. 3. The magnetization texture of a liquid ferromagnetic thin film elongating under applied field. The magnetization 
texture in zero field is going counterclockwise. Values for dimensionless field h* (see Eq. (|lf])) are 0, 0.42 and 0.68. 



9 




20 40 60 80 100 120 140 



2r /A 

FIG. 4. Plots of e/A vs 2ro/A for ferromagnetic and paramagnetic droplets for a fixed value of Ho = 20 G. 
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